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When a spin-polarized current flows through a ferromagnet, the local magnetization receives a 
spin torque. Two consequences of this spin torque are studied. First, the uniformly magnetized 
ferromagnet becomes unstable if a sufficiently large current is applied. The characteristics of the 
instability include spin wave generation and magnetization chaos. Second, the spin torque has 
profound effects on the structure and dynamics of the magnetic domain wall. A detail analysis on 
the domain wall mass, kinetic energy and wall depinning threshold is given. 
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I. INTRODUCTION 

There are many interesting phenomena generated by 
electronic current. In this paper, we focus our discussions 
on two effects: domain wall dynamics and magnetization 
instability. Recently, both topics have received consider- 
ably interests in experiments GIB El, El E l EI B E I 0, 0, 
EH and in theories [H [H MJUl HIllE Ellll It is 
shown that the current is able to displace magnetic do- 
main walls in a spin valve 0, in a constricted nanowire 
0| , in U-shaped (or L-shaped) nanowires @> , in a ring 
structure Q and in zigzag wires [j| . The current also gen- 
erates dynamic domain walls (time-dependent wall struc- 
ture) via the generation of spin waves in the uniformly 
magnetized ferromagnet. The point-contact experiments 
on a single layer ferromagnet p,ITol[Tll| are likely associ- 
ated with the interface spin waves excitations [12 ]. 

The physics of these experimental results is believed to 
be the consequence of spin angular momentum transfer, 
originally proposed by Slonczewski and Berger |2^| 
in magnetic multilayers. The extension of Slonczewski's 
model to a homogeneous ferroma gne t has been carried 
out by several groups [iH ITU Il6l 1l7j. Recently, we in- 
troduce two spin torques in ferromagnets [Tsj . With our 
spin torque model, most of the experimental observations 
can be quantitatively analyzed. This paper is organized 
as follows. In Sec. II, we briefly review the mechanism 
of the spin torque in a single ferromagnet. In Sec. Ill, 
we study the bulk and surface spin waves by applying 
the spin torque model. Several key characteristics of the 
magnetization instability driven by spin torques are pre- 
dicted. In Sec. IV, we analyze the domain wall dynamics 
with and without defects. We also construct an ana- 
lytic 1-D model to help understanding the key features 
in simulations. Finally, an outlook for the effects of spin 
torques in ferromagnets is presented in Sec. V. 

II. MECHANISM OF SPIN TORQUES IN 
FERROMAGNETS 

The spin-polarized current does not interact with spa- 
tially uniform magnetization, except the classical Zeeman 
coupling between the current-induced magnetic field and 



the magnetization. In a real ferromagnet, the magnetiza- 
tion is not spatially uniform due to the presence of vari- 
ous magnetic interactions, which tend to break the ferro- 
magnet into domains. Even in a small particle, where 
the single domain assumption is approximately valid, 
the thermal fluctuation leads to a time-dependent non- 
uniform magnetization. In response to the spatially non- 
uniform magnetization, the spin current will be position- 
dependent and the spatial varying spin current produces 
the mechanism of the spin angular momentum transfer. 
To determine the spin transfer torque, one requires si- 
multaneously solving the magnetization dynamics and 
the time-dependent spin current, which would be a very 
complex problem. Fortunately, the dynamics of the mag- 
netization is much slower than that of the transporting 
electrons, then it is a good approximation to work on the 
transport equation with magnetization " frozen" spatially 
and temporally. Similar to the calculation of electron dy- 
namics in an atom: the electronic structure is calculated 
by freezing the dynamics of nuclei. 

Most of the theories of spin torques are built on 
the above approximation. Additional assumptions are 
needed to derive analytical expressions of the spin torque. 
Berger 23] first investigated the current-induced domain 
wall motion and introduced the "domain drag force" by 
using an intuitively physics picture that the current can 
drag the domain wall moving along the path of the cur- 
rent flow via an s-d exchange interaction. Since the spin 
torque was not mathematically formulated in this work, 
it is unclear how to solve the problem of domain wall 
dynamics for a realistic magnetic wire by using this intu- 
itive approach. Bazaliy et al 0] proposed a spin torque 
model in a ferromagnet within the ballistic transport 
model for half-metallic materials and they found that 
the spin torque is t a (j e ■ VM) where j e is the elec- 
tric current. The essential assumption is that the spin 
polarization of the current is parallel to the local mag- 
netization, i.e., an adiabatic approximation. The above 
explicit expression of the spin torque can be immediately 
combined with the well-known LLG equation to calculate 
the response of magnetization to the spin torque. Wain- 
tal and Viret |16T ] has extended this approximation by 
relaxing the adiabatic approximation. They have shown 
that the spin polarization of the current is not parallel 
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to the local magnetization in a ballistic transport model. 
However, no clear mathematical expression is given to 
implement this additional spin torque to the LLG equa- 
tion. Tatara and Kohno |17| proposed two spin current 
effects: an adiabatic torque mentioned above and a mo- 
mentum transfer torque. The momentum transfer effect 
originates from the momentum scattering by a domain 
wall and it is proportional to wall resistivity. This mo- 
mentum transfer effect is negligible except for very thin 
walls. 

We have proposed the spin torque by evaluating the re- 
sponse of the conduction electron spins in a spatially and 
temporally varying magnetization M(r, t) in the semi- 
classical transport theory ^| . The spin torque has been 
formulated in the following form 



bj(j e ■ V)M - -g-M x (j e 



V)M 



(1) 



where j e is the unit vector in the direction of the cur- 
rent flow, bj — Pj e fiB I 'eM s , P is the spin polarization 
of the current, M s is the saturation magnetization, /is is 
Bohr magneton, cj — and £ is a dimensionless con- 
stant that describing the degree of the nonadiabaticity 
between the spin of conduction electrons and the local 
magnetization. For a typical ferromagnet (Ni,Co,Fe and 
their alloys), £ is within a range of 0.001 ~ 0.05. The 
physical interpretation is that the bj term is an adia- 
batic spin torque, describing the adiabatic process of the 
non-equilibrium conduction electrons, i.e., the direction 
of the spin polarization of current is parallel to the lo- 
cal magnetization within a domain wall; the cj term is a 
non-adiabatic torque, which is related to the spatial mis- 
tracking of spins between conduction electrons and local 
magnetization. 

The advantage of expressing the spin torques in the 
form Eq. (1) is that one can readily generalize LLG equa- 
tion in the presence of the spin current, 



— = - 7 M x H e// 



Q' <9M 
— M x — 

M s dt 



-*t+t<> ( 2 ) 



where 7 is the gyromagnetic ratio, and H e ff is the ef- 
fective magnetic field including the external field, the 
anisotropy field, magnetostatic field, and the exchange 
field, and a is the Gilbert damping parameter. The 
current-driven magnetization dynamics will be studied 
by solving this generalized LLG equation in various sit- 
uations. 



III. SPIN WAVE EXCITATIONS 



where Hk is the anisotropy field, alone the x-axis, 
A is the exchange constant, and we include a self- 
demagnetization field 4irM z of the film. The initial mag- 
netization saturates in the direction of the magnetic field 
H e along the x-axis in the plane of the layer. 

When one applies a sufficiently large current along the 
x-axis, the uniformly magnetized film becomes unstable. 
To see this, we consider a small deviation of the magne- 
tization vector from the easy axis e Xl 



M = ALe, 



= i(a)t+k-r) 



(4) 



where <5m is a small vector. Inserting Eqs. (3) and (4) 
into Eq. (2), and keeping only the terms linear in <5m, 
we obtain two linearized equations for Sm y and Sm z . A 
secular equation is then established for the spin wave 
frequency u> and the spin wavevector k. 



- bjk x ) 
A 2 



- bjk x ) 



Sm y 
Sm z 



(5) 



where we have defined A\ 



-7 I W 



4ttM, + H 



icjk x and A2 = 7 ( jr-k 2 + H ) + iau — icjk 



which k 2 = k 2 + k 2 + k 2 and H = H e + H K . 

x y z 

The above linearized equations have a non-zero solu- 
tion if and only if 



det 



-i((j-bjk x ) A 1 

A 2 -i(u> - bjk x ) 







(G) 



The above equation establishes the relation between the 
spin wave frequency u> and the wavevector k. For a real 
wavevector k, u is a complex number. If the imagi- 
nary part of to becomes negative, Eq. (4) is exponentially 
growing with time. In this case, the uniform magneti- 
zation becomes unstable. One may define the critical 
current such that Im cj = 0. From Eq. (6), we find 
Re uj = cjk x /a and 



cj/a 



k T \ \ M. / V M, 



(7) 



Equations (7) gives out the current density required to 
generate spin waves with a given wave vector k. Inter- 
estingly, the minimum current density does not occur at 
the uniform mode of k = as in the ordinary spin wave 
excitations, rather the spin wave with a finite wavevector 
is first excited by the current. To see this, we minimize 
bj — cj/a in Eq. (7) with respect to the wavevector, and 
we find the minimum current density occurs at 



In this section, we apply our generalized LLG equation, 
Eq. (2) to study the spin wave excitations in a single layer 
magnetic film. The local effective field is 



A, .... 



M. 



-^^{AnM s + H)H. 



2A 



(8) 



and k y = k z = 0. The corresponding wavelength is, 



H K M X 

f-pff — TT e i 



2A 

mF 



V 2 M - 4TrM z e z + H e e x (3) 



H)H] 



-1/4 



(9) 
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The new length scale given above places a severe limi- 
tation on micromagnetics: one needs to choose mesh size 
to be smaller than A c when the current density exceeds 
the critical current density, in order to correctly capture 
the current-driven effects in the simulation. For certain 
parameters, particularly at large magnetic fields, A c can 
be comparable or smaller than the exchange length. 

Inserting Eq. (8) back to Eq. (7), we find that the 
minimum current density for the instability is, 



c,i/a\ mm =lJjf + VH + 4ttM s J . (10) 



Note that we have assumed a spin wave spectrum in 
the form of Eq. (8) which represents a spin wave in an 
isotropic infinite medium. If we take £ = a the spin wave 
instability does not occur for any large spin currents. 

In experiments, the current may not uniform across 
the sample. For example, the current density is very 
large at the contact area in the point-contact experiments 
and it becomes negligibly small away from the contact 
area. In this case, the excited spin waves are confined 
in a narrow region near the interface. To access the spin 
wave instability in this situation, we consider a surface 
spin wave mode 



M = M s e x + 5 me - x / K e l( - UJt+k y y+k * z \ 



(11) 



where k is the penetration length. 

By inserting it into Eq. (2) and by repeating the deriva- 
tion similarly, we find the instability Im lo = occurs at 



-(bj+acj) = jan (j^-(k 2 y + k 2 ) + H + . (12) 

Note a negative sign on the left-hand side; it indicates 
that the surface spin waves can only be generated by the 
current flowing along one direction. This result is consis- 
tent with the point-contacted experiments [lfj , in which 
the current flowing from the magnetic layer to the non- 
magnetic tip is able to excite spin waves. The magnitude 
of the critical current for surface spin wave is propor- 
tional to the damping parameter a and the penetration 
length k. We note that the exchange field at the surface 
for the localized spin wave is (2Aao/M|)V x M where ao 
is the lattice constant so that the exchange field inside 
the layer {2A/M 2 )V 2 X M will be compensated. This is 
the reason that the right hand side of Eq. (12) does not 
contain the k 2 term. 

An estimation on the magnitude of spin wave instabil- 
ity can be readily done by using the materials parameters 
of Co: 7 = 1.9 X 10 7 (Oe)- 1 s- 1 , AttM s = 1.8 x 10 4 Oe, 
H K = 500 Oe, M s = 14.46 x 10 5 A/m, A = 2.0 x 
10~ n J/m and P = 0.35. If one takes the damping 
parameter to be 0.01, the non-adiabaticity ( — 0.02, the 
penetration length of 5 nm, and k y = k z « for the 
long wave length limit, we find that the critical current 
at H e = Oe is = 1.12 x 10 10 A/cm 2 for the bulk 



spin wave and J^, ° e = 6.4 x 10 7 A/cm 2 for the sur- 
face mode, which is two orders of magnitude smaller than 
that of bulk spin wave. 

The fact that the large enough current density can ex- 
cite spin waves with different wave lengths raises a ques- 
tion: what is the magnetization state at the large cur- 
rent density? Shibata et al |19| predicted that a static 
multi-domain state can be formed if the current exceeds 
a second critical current higher than the critical current 
density defined above. By using Eq. (16), we find that 
the static domain configuration is unlikely to form. This 
is because the domain will be moving with an average 
velocity determined by the non-adiabatic torque cj/a 
[iH [20j seen in next section. We had reported earlier 
that a large current density drives the uniform magneti- 
zation into spatially and temporally chaotic motion |2l| . 
The detail analysis on this transition will be given else- 
where. 



IV. DOMAIN WALL MOTION 

One of the most interesting predictions on the domain 
wall motion is that a steady domain wall velocity is inde- 
pendent of bj. Rather the steady state velocity is solely 
determined by the non-adiabatic torque and the external 
magnetic field. We show the proof below. 

In a steady state motion, one can assume the magne- 
tization vector M = M(x — v x t,y,z). Then, dM/dt = 
— v x dM/dx. For a current applied in x-axis, the LLG 
equation, Eq. (2), is thus 



{bj + v x ) 



dM 

dx 



av x 



AL 



ox 



7 M x H e// (13) 



By performing inner-product of M x dM/dx with the 
above equation, one has 



av x 



M x 



dM 

dx 



= 7 H e// 



dM 

dx 



(14) 



We now integrate the above equation over the entire do- 
main wall. We find 

2 



e.j 



dV 



dM 



dx 



dx v x dt 



EdV 



(15) 

where E is the energy density that is defined dE = 
—H e // • dM, and we have also used the fact that M 
and dM/dx is perpendicular. The right-hand side of 
the equation becomes ±2SoM s H e because the rate of 
the total energy change for a uniformly moving wall is 
2SoM s H e v Xl where So is the cross section perpendicular 
to the moving wall. Therefore, the velocity is 



v, 



jH e W 



a a 

where we have defined the wall width W 



w- 1 = 



1 



2S M 2 



dV 



dM 



dx 



(16) 



(17) 
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The striking conclusion from Eq. (16) is that the wall 
velocity in the absence of the magnetic field is exactly 
—c,j/a for any type of walls as long as the wall is moving 
at a constant velocity. Without the curr ent. Eq. (16) be- 
comes the well-known Walker's velocity |24|. We should 
point out that the wall width, Eq. (17), is weakly de- 
pendent on the magnetic field and the adiabatic torque 
bj 0], an d thus the terminal wall velocity will slightly 
depend on bj in the presence of the field. 

The wall velocity given by Eq. (16) breaks down when 
the wall is not moving with a constant velocity. To de- 
scribe a non-uniform wall motion, it is necessary to pos- 
tulate an approximate wall structure. Here, we follow 
Walker's procedure [24| by introducing two polar angles 
8 and <p in the form mimic a transverse domain wall [25l | . 
By placing the explicit form of 9 into Eq. (2), we find 

^ = 7 (if e - 47raAf s sin p cos tp) + abj ^ ° J (18) 

av x (t) + W^ = 1 H E W-c J (19) 
at 

where ip is the out-of-plane angle of the magnetization 
vector, v x (t) is the velocity, and we have discarded a 2 
terms. 

To further gain insights of the wall velocity, we as- 
sume that the distortion of the wall is small during the 
wall motion. Then, we treat the wall width W as time- 
independent constant, and replace sin t/j cost/? by ip in 
Eq. (18). By differentiating Eqs. (18) and (19) with re- 
spect to t, and by eliminating ip in the resulting equa- 
tions, we find the equation for the wall velocity 

± dv x (t) 2M s a , s 2M S , s 

dt 

where we have defined the domain wall effective mass per 
unit area as j2(J , 



2ij,\-1 



m* = (2-kYW) 



(21) 



Equation (20) reveals that the domain wall can be treated 
as a particle with an effective mass to*, subjecting to a 
friction force (second term) and external forces from the 
magnetic field (third term) and the current (fourth term) . 
Interestingly, the adiabatic torque does not contribute to 
the current-driving force. The solution of Eq. (20) is 



Vx(t) = VCe 1 



a 



(22) 



where t = (4tt Mg'ya) is the relaxation time and C is 
a constant determined by the initial condition. If one 
assumes that the field or the current is applied at t = 0, 
i.e., <p(Q) = 0, we find, from Eq. (18) and (19), v x (0) = 
—bj. By using this initial value of the velocity, Eq. (22) 
becomes 



v x (t) 



cj - jWH e 



a 



(l - e-*/ T ) - bje-^ (23) 



One may understand the above two terms as follows. The 
first is a "translational" velocity and one might introduce 
a kinetic energy 



1 , 2 (CJ - jWH e ) 2 

E = -TO v x (t) = — - — — 

2 x w 4ir"f 2 Wa 2 



1-e 



-t/T 



(24) 



The time scale to establish a fully accelerated motion is 
determined by r. For a Co wire, t = 0.3 ns if we take 
a = 0.01. The second term of the Eq. (23) is related 
to the domain wall distortion, and therefore, it does not 
contribute the mobility of the domain wall. The displace- 
ment of the domain wall is x c — J Q v x (t)dt. 

We illustrate below that the concept of the "trans- 
lational" kinetic energy is helpful in understanding the 
domain wall motion in a nanowire containing defects or 
pinning centers. For example, let us consider a domain 
wall propagation from one place (A) to another place 
(B), shown in Fig. 1(a). Between A and B regions, there 
is a distribution of defects. We show that the minimum 
current-density or the minimum magnetic field required 
to overcome these pinning centers depends on the relative 
distance between the defects and the domain wall. If a 
strong defect is located far away from the initial domain 
wall position (A), the domain wall is able to develop its 
fully kinetic energy, and the domain wall is able to pass 
through the pinning center when the kinetic energy is 
larger than pinning potential. On the other hand, if the 
defect locates near the initial wall center, the kinetic en- 
ergy is small when the domain wall encounters the defect 
and thus a relative weak defect potential is able to trap 
the domain wall. To be more quantitative, we perform 
the following calculation. 

We model a defect by introducing an artificial lo- 
cal anisotropy pinning source in an otherwise perfect 
nanowire. The anisotropy of the defect is Hd = 4Hk 
which is 4 times larger than the anisotropy of the wire, 
Hk = 500 Oe. We choose a one-dimensional model, i.e., 
the magnetization does not vary in the direction of width 
and thickness, to illustrate our points. The easy axis of 
the wire and the defect, and the current or magnetic field, 
are all along x-axis. The mesh size in x-direction is 4 nm. 
Before one turns on the current, a stationary Neel wall is 
centered at x = 0, and the defect is located Xd from the 
center of the wall. At t — 0, a current is turned on and 
we calculate the domain wall motion afterwards. 

First, we have placed the defect far away from the re- 
gion A, i.e., Xd 3> J v x (t)dt, we find that the critical 
current is proportional to the damping parameter in the 
absence of magnetic field, consistent with Eq. (24). When 
we vary the exchange constant which effectively changes 
the domain wall width, we find that the critical current 
scales as VW, again, consistent with Eq. (24). We now 
vary the defect position. In Fig. 1(b), we show the crit- 
ical current as a function of the defect position. As ex- 
pected, the critical current decreases as one places the 
defect away from the original wall. The initial increase 
of the critical current at small distance comes from the 
adiabatic torque, bj. When the defect is very near the 
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FIG. 1: (a) Schematic pinning source in a Co nanowire, (b) 
The minimum spin current required for a wall propagation 
through a pinning as a function of pinning position Xd- The 
parameters of Co are 4-kM s = 1.8 x 10 4 Oe, H K = 500 Oe, 
Ms = 14.46 x 10 5 A/m, A = 2.0 x 10 -11 J/m, H e = 0, 
C = 0.02 and a = 0.006 



original wall, the displacement of the wall by the bj term 
exceeds Xd and thus the wall overcomes the defect po- 
tential by the adiabatic torque. Therefore, at the small 
distance, the adiabatic torque is more effective than the 
kinetic energy involved. 



V. OUTLOOK 

There are many fundamental issues on current-induced 
effects. We intend to list a few of them below. 

On the fundamental mechanism, the spin torque given 
by Eq. (1) is valid up to the first derivative in space. 



When the domain wall thickness is extremely small in 
some nanoconstriction, the high order terms become im- 
portant. A better quantum mechanical treatment for the 
spin torque is required to re-formulate the spin torque in 
this case. 

In a large field or a large current, the wall is usually 
moving quite irregularly. Non-uniform wall motion ap- 
pears: the wall can be bounced back and forth in trans- 
verse directions, and the wall changes its structure. For 
example, the current can drive vortex wall motion rela- 
tively easier than the transverse wall in the presence of 
defects. However, the vortex wall is usually unstable dur- 
ing its motion and the vortex wall tends to change into 
a transverse wall. These detailed complications have to 
be studied via numerical solutions of the LLG equation 
[20ll27| . To model a realistic device, an extensive numer- 
ical effort is required. 

The effect of Joule heating by the current has not been 
quantified. The experimental results show that the sam- 
ple temperature increases less than 5 K when the cur- 
rent density is order of 10 5 ~ 10 6 A /cm 2 3], while the 
temperature is dramatically increased when the current 
density is 10 7 A/ cm 2 plf. These studies are for a steady 
current density. In most of the devices, it is much desir- 
able to use a pulsed current source. A thorough study on 
the Joule heating on the amplitude and duration of the 
current seems very important. Frequently, the observed 
critical current is smaller than the optimal theoretical 
value by a factor of two. Whether it is due to heating is 
unclear at the present time. 

The finite temperature theory on the domain wall mo- 
tion is lacking. It would be interesting to study the life- 
time of a trapped domain wall and it takes for the cur- 
rent to depinning the wall at finite temperature. Similar 
to the effect of the finite temperature in current-driven 
switching in spin valves, we expect that the role of tem- 
perature has profound effects on domain wall motion. 

The research was supported by NSF (ECS-0223568 and 
DMR-0314456). 
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